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ABSTRACT 

We studied Faraday rotation measure (RM) in turbulent media with the rms 
Mach number of unity, using isothermal, magnetohydrodynamic turbulence sim- 
ulations. Four cases with different values of initial plasma beta were considered. 
Our main findings are as follows. (1) There is no strong correlation between the 
fluctuations of magnetic field strength and gas density. So the magnetic field 
strength estimated with RM/DM (DM is the dispersion measure) correctly rep- 
resents the true mean strength of the magnetic field along the line of sight. (2) 
The frequency distribution of RMs is well fitted to the Gaussian. In addition, 
there is a good relation between the width of the distribution of RM/RM (RM 
is the average value of RMs) and the strength of the regular field along the line 
of sight; the width is narrower, if the field strength is stronger. We discussed the 
implications of our findings in the warm ionized medium where the Mach number 
of turbulent motions is around unity. 
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Introduction 



The warm ionized medium (WIM) is one of the major gas components of our Galaxy. It 
is a diffuse ionized gas with temperature T ~ 8000 K, scale height H ^ 1 kpc, and average 
density n ~ 0.03 cm~^, and occupies approxima tely 20% of the volume of the disk in the 



Galaxy (e.g., Reynolds 1991 : Haffner et al. 19991). The measur ed width of the Ha line from 



the WIM is in the range of 15 to 50 km s~^ (ITufte et al.lll999l ). Observations indicate that 
the WIM is in a turbulent state (see below). Assuming that the typical width of Ha line and 
temperature are 30 km s~^ and 8000 K, respectively, the non-thermal, turbulen t velocity 
would be about 13 km s~^, as is calculated, for instance, with Equation (1) of [Reynolds 
( 119851 ). So the sonic Mach number of turbulent motions, Mg, would be ar ound unity in the 



WIM. It is smaller than Mg of the cold neutral medium, which is a few (e.g.. lHeiles &: Troland 



20031 ). and Mg of molecular clouds, which is > 10 (e.g.. lLarsonlll98ll ) 



The best evidence for turbulence in the inters tellar medium (ISM) w hich includes the 
WIM comes from the power spectrum presented in [Armstrong et al.l (119951 ). It is a compos- 
ite power spectrum of electron density collected from observations of velocity fluctuations 
of the interstellar gas, rotation measures (RM), dispersion measures (DM), interstellar scin- 
tillations, and others. The spectrum covers the spatial range of ~ 10 



10 



10 



20 



cm, 



and 



remarkably the whole spectrum is fitted to the power spectrum of Kolmogorov turbulence 
with slope —5/3. In additi on, it has been recently reported that the Ha emission measure for 
the WIM ([Hill et al.l 120081) and the densities of the diffuse ionized gas and the diffuse atomic 
gas (iBerkhuijsen fc Fletcherl l2008l ) follow the lognormal distribution. The lognormality in 
density (or colum n density) distributions can be regarded as another signature of turbulence 



in the WIM (e.g.. lVazquez-Semadenilll994j : lElmegreen fc Scaldl2004l . and references therein). 



The interstellar magnetic field that is pervasive in the Galaxy plays an important role 
in the dynamics of the ISM, star formation, acceleration of cosmic rays, etc. It has the 
energy density c omparable to those of turbulence and cosmic rays, as well as the thermal 



energy density (jSpitzerl Il978l ). The information on the field has been obtained through 
observations of Zeeman splitting, polarized thermal emission from dust, optical starlight 
pola rization, radio synchrotron emission, and the Faraday rotation of polarized radio sources 



see. 



Han &: Wielebinskil 120021 . for a review). Of them, the last method can be used to study 



the magnetic field in ionized media such as the WIM. It estimates the mean strength of the 
magnetic field along the line of sight, weighted with the electron density rig, 

Jn^B^^ds _ RM 
= "/MT = DM' 

where B \ \ = B cos 6 and 9 is the angle between B and the line of sight. Several authors (e.g. , 
Han et al.l[l998l : llndrani fc Deshpand^ Il999l : IPrick et al.lboOll : iHan et al.l[2006l : ISecklboOTl ) 
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have used the RMs and DMs of pulsars to reproduce the large-scale magnetic field of our 
Galaxy. 

Recently, the distributions of RMs along many contiguous lines of sight have been 
obtained in multi-frequency polarimetric observations of the diff use Galactic synchrotron 
background (iHaverkorn et aL l2003l . l2004j : ISchnitzeler et all boooh and the Perseus cluster 



(Ide Bruyn fc Brentjend l2005l ). While the peak in the frequency distribution of the RMs 
refiects the regular component of the magnetic field, the spread should refiect the turbulent 
component. This means that if the distribution of RMs is observed, the spread provides 
another way to quantify the magnetic field in turbulent ionized media such as the WIM. 

Motivated by the importance of RM in the study of the magnetic field in the WIM, in 
this paper we study RM in turbulent media with the rms Mach number of unity. Simulations 
are outlined in Section 2. Results are presented in Sections 3 and 4. 



Simulations 



We performed three-dimensional simulations of isothermal, compressible, magnetohy- 
drodynamic (MH D) turbulence, using a code based on the total variation diminishing scheme 
(IKim et al.lll999l ). The fact that the te mperature of the WIM is in a relatively narrow range 
of 6000 - 10000 K (IHaffner et al.lll999l ) would justify the assumption of isothermality. There 
are two parameters in our simulations: the root-mean-square (rms) sonic Mach number, M^, 
and the initial plasma beta value, jSo- Hereafter, the subscript "0" is used to denote the initial 
values. The first parameter indicates the level of turbulence. In this paper, we focus on the 
turbulence with the rms Mach number of unity, and so set Mg = 1. The second parameter 
tells the strength of the initially uniform magnetic field, or the regular field, Bq. In order to 
explore the effect of regular field, we varied the value of jSo from 0.1, 1, 10, to 100. If we take 
8000 K and 0.03 cm~^ as the representative values of the temperature and electron density, 
assuming that hydrogen is completely ionized, helium is neutral, and the number ratio of 
hydrogen to helium is 0.1, we have Bq = 1.3(l//3o)^/^(T/8000 Ky/'^{njOm cm-3)V2 ^g. 
So the initial magnetic field strength in our simulations corresponds to 4.1, 1.3, 0.41, 0.13 
jj,G for /?o = 0.1, 1, 10, and 100, respectively. 

Simulations were started with Bq along the x-direction in a uniform medium of density, 
Pq. The grid of 512^ zones w as used for the peri odic computati o nal b ox of size L. Turbulence 
was driven with the recipe in lStone et al.l (Il999l ) and lMac Lowl (Il999l ): velocity perturbations 
were drawn from a Gaussian random field determined by the top-hat power distribution in 
a narrow wavenumber range of {2n/L) < k < 2{2tt/L), and added at every At = O.OOlL/c^. 
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The amplitude of the perturbations was tuned in such a way that Mg = frms/cs became 
nearly unity at saturation. Here, is the isothermal sound speed and Vj-ms is the rms 
velocity of the resultant turbulent flow. In simulations frms initially increased and became 
saturated at t ~ {l/2)L/cs, and we ran the simulations up to t = 2L/cs, two sound crossing 
times. At saturation, the average strength of magnetic field reached 4.2, 1.45, 0.73, and 0.51 
/iG for Po = 0.1, 1, 10, and 100, respectively. The amplification factor was larger for the 
simulations with larger Pq or weaker Bq. 



3. Correlation between B and p 



We first check whether {B\\) in Equation (1) correctly reproduces the unbiased strength 
of the magnetic field along the line of sight. That is true, only if there is no correlation 
between B and He (or the gas density p). If the correlation is positive, {B\\) in Equation (1) 
overestimates the magnetic field strength, while if negative, it underestimates the strength. 



Ostriker et al 



2005; Mac Low et al. 



2005 



20011: 



It is known fro m observations (ICrutchei 
cal sim ulations (e.g., 



199S 



Padoan &: Nordlund 1999) and numeri- 



Passot fc Vazquez-Semadeni 



2003 



Balsara fc Kim 



Burkhart et al. .200ft) that in the highly compressible, super- 



sonic, molecular cloud environment, the mag netic pressure (or B) a nd the gas pressure (or 
p in isothermal gas) are positively correlated. iBurkhart et al.l (120091 ). on the other hand, re- 
ported a weak negative correlation for Mg = 0.7 and (3q = 2 (their model 1). Figure 1 shows 
the correlation between B and p in our simulations. The nearly circular shape of contour 
lines indicates weak correlation for Mg = 1. To quantify it, we calculated the correlation 
coefficient 



r{B,p) 



(2) 



where B and p are the average values of B and p. The values of r at the end of simulations 
were 0.01, -0.12, -0.06, and 0.16 for /9o=0.1, 1, 10, and 100, respectively. Relatively small 
values of r confirm that the correlation is weak for Mg = 1, regardless of (3o- 

Weak correlation means that the RM field in Equation (1) should correctly represent 
the true magnetic field. To test it, we compare the strengths of two fields in Figure 2; the 
frequency distributions of the true mean strength of the magnetic field along the line of sight 
(presented with solid lines) and the magnetic field strength calculated with Equation (1) 
(presented with dotted lines) coincide quite well. So we argue that the systematic bias, due 
to a correlation between B and ng, in the estimation of magnetic field strength with RM 
would be insignificant in the turbulent media with Mg = 1. 
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Beck et al.l ( 120031 ) pointed that the discrepancy in the strengths of the regular Galactic 



magnetic field estimated with RM and synchrotron emissivity could be reconciled, if the cor- 
relation between B and rie is negative in the WIM and so the RM field was underestimated. 



Beck et al.l ( 120031 ) postulated the pressure equilibrium, which results in the negative corre- 



lation between B and p. If the pressure equilibrium is maintained, the combined pressure 
of gas and magnetic field, Ptot = -Pgas + -Pmag, should exhibit a narrow distribution. Figure 
3 shows the frequency distribution of Ptot- Our simulations show broad distribut ions rather 



than p eaked distributions. That is, our results do not support the prediction of iBeck et al. 
(120031 ). However, we caution that we considered only the turbulent media with the rms 
Mach number Mg = 1. On the other hand, the WIM may have Mg which is not exactly 
unity. So we nee d to further investigate turbulent media with Mg ^ 1, before we exclude 
the postulation of iBeck et al.l (120031 ). We leave it as a future work. 



4. Frequency Distribution of RMs 

We then look at the frequency distribution of RMs. Figure 4 shows the probability 
distribution of RM/RM for different jS^s as well as for different ^'s. Here, RM is the average 
value of RMs. (The distribution was also calculated for 6 = 83°, but not shown in the 
figure for clarity. The distribution for 6 = 83° is used in Figure 5.) The fit to the Gaussian 
is also shown. A noticeable point in the figure is that the distribution of RM/RM is very 
well fitted to the Gaussian. The go odness-of-fit is be t ween 0.89 and 0.99 0. This result 



is consistent with the observations of iHaverkorn et al.l (120031 . |200J). They took the multi 



frequency polarimetric images of diffuse radio synchrotron backgrounds in the constellation 
Auriga and Horologium, and obtained RM maps. The distribution of the observed RMs is 
also fitted to the Gaussian. 

Another noticeable point in Figure 4 is that the distribution is more widely spread for 
larger j3o and larger 6. It hints at a possible correlation between the spread in the distribution 
of RM/RM with the strength of the magnetic field along the line of sight. Figure 5 shows the 
relation between the full width at half maximum, VFfwhm, of the distribution of RM/RM 
shown in Figure 4 and the strength of the regular field along the line of sight. Poll = Bq cos 6. 



^We used the coefficient of determination, i?^, as an indicator of goodness-of-fit, wliere i?^ = 1 if the fit 
is perfect and i?^ < f if the fit is less perfect. 
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The best fit for the relation B is 



Bow = (2-45 ± 0.3) X W^f"™"'' f^G. 



(3) 



Note that the relation is for T = 8000 K and ne = 0.03 cm~^ as the representative values 
of the temperature and electron density, and scales as (T/8000 K)^/^(r;,e/0.03 cm~'^)^/^ for 
other values of T and rig. The empirical relation in (3) may provide a handy way to quantify 
the strength of the regular field along the line of sight in regions where the map of RMs has 
been obtained along many contiguous lines of sight with the multi-frequency polarimetric 
observations (see Introduction). The accompanying map of DMs is not necessary for it. 

The estimation of magnetic field strength with the relation in (3), however, should be 
done with caution because of its limitations: First, the relation is valid only for Mg = 1. 
Second, no stratification of gas and magnetic field was assumed, which is apparent in the 
polarization maps of a larger area of our Galaxy. Third, no source emitting polarized light 
in the medium between background polarized continua and us was considered. 

From a physical point of view, the broadening of the width of RM distribution is due to 
fiuctuating magnetized gas. We can easily expect that, for a given regular field, the width 
would be larger if the rms Mach number, Mg, is larger. So Wfwhm should be a function of 
not only _Bo|| but also Mg. In order to quantify the relation of VFfwhm versus Bq\\ and Mg, we 
need simulations that cover the parameter space of /5o and Mg. Low-resolution simulations 
with Mg = 0.5 and 2.0 show that the probability distribution of RM/RM is still fitted to the 
Gaussian, and Wfwhm can differ by a factor of two to three if Mg differs by a factor of two. 
We leave the report of the Mg ^ 1 results, including the correlation between B and p, from 
high-resolution simulations future work. 

Nevertheless, we can try to apply the relation in (3) to RM observations in the WIM, 
assuming Mg = 1 there. For i nstance, from the f igures of the frequency distribution of RMs 
in Auriga and Horologium in iHaverkorn et al.l (120031, I2004J ) , we read Wfwhm — 2 and 7, 
and get Bq\\ ^ 0.6 and 0.1 AtG for Au riga and Horologium, respectively (if Ue = 0.016 cm~'^ 



is used as in IHaverkorn et al. J 200 



values estimated in Haverkorn et al 



[) ). Th ese are close to, but somewhat larger than, the 
(120041 ) ■ ~ 0.42 and 0.085 /iG, respectively. However, by 
considering the limitations of the relation in (3) and the uncertainty in the values of Wfwhm 
which we read from figures, we regard the agreement to be fair. 



Currently the number of the synchrotron backgrounds, which can be used for the obser- 



^We note that in the fit, there is a trend that the blue symbols for /3o = 0.1 are mostly above the fitting 
line, while the black symbols for po = 100 are mostly below the line. This tells that Pq may enter the relation 
as a secondary, but less prominent, parameter. 
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vation of RM maps, is still limited due to the relatively low sensitivity of present-day radio 
telescopes. However, the new-generation radio telescopes, such as LOFAR (Low Frequency 
Array) and SKA (Square Kilometer Array) with much higher sensitivity, will certainly pro- 
vide RM maps in much larger portions of the sky. Then, the relation like the one in (3) will 
become a useful diagnosis. 
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Fig. 1. — Contour plots showing the correlation between the magnetic field strength and the 
gas density at the end of fours simulations. The contour levels are from 10% to 90% of the 
peak value with the interval of 10%. 
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Fig. 2. — Frequency distribution of (a) {B^) and (b) {By), normalized with the initial 
magnetic field strength, Bq. The solid lines present the true mean strength of the magnetic 
field along the line of sight, and the dotted lines present the magnetic field strength calculated 
with RM/DM, along 512^ lines aligned the (a) x and (b) y-axes. The blue, red, green, and 
black lines are for the simulations with (3q = 0.1, 1, 10, and 100, respectively. 
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Fig. 3. — Frequency distribution of the total (gas and magnetic) pressure. The blue, red, 
green, and black lines are for the simulations with jSo = 0.1, 1, 10, and 100, respectively. 
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RM/RM 

Fig. 4. — Probability distribution of RMs normalized with the average value. The solid lines 
are the distribution, and the corresponding dashed lines are the Gaussian fit. The blue, red, 
green, and black lines are for (3q = 0.1, 1, 10, and 100, respectively. Different panels are 
for different values of 6 (the angle between the direction of the regular field and the line of 
sight). 
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Fig. 5. — Relation between the full width at half maximum of the distribution of RM/RM 
(W^fwhm) and the strength of the regular field along the line of sight (-Bo||)- Blue, red, green, 
and black symbols are from the simulations with = 0.1, 1, 10, and 100, respectively (for 
five different ^'s including 6 = 83°). The error bar shows the error in the Gaussian fit in 
Figure 4. The solid line shows the best fit of the relation between VFfwhm and i?o||- 



